Under the influence of alcohol: The effect of ethanol and methanol on lipid 
bilayers 

Michael Patra 

Biophysics and Statistical Meclianics Group, Laboratory of Computational Engineering, Helsinki University of Technology, 
P. O. Box 9203, FIN-02015 HUT, Finland 

Emppu Salonen, Emma Terama, and Ilpo Vattulainen 

Laboratory of Physics and Helsinki Institute of Physics, Helsinki University of Technology, P. O. Box 1100, FIN-02015 HUT, 
Finland 

Roland Faller and Bryan W. Lee 

Department of Chemical Engineering and Materials Science, University of California-Davis, One Shields Ave, Davis, 
CA 95616, USA 

Juha Holopainen 

Department of Ophthalmology, University of Helsinki, Finland and Helsinki Biophysics & Biomembrane Group, 
Institute of Biomedicine, University of Helsinki, Finland 

Mikko Karttunen 

Biophysics and Statistical Mechanics Group, Laboratory of Computational Engineering, Helsinki University of Technology, 
P. 0. Box 9203, FIN-02015 HUT, Finland 

Extensive microscopic molecular dynamics simulations have been performed to study the effects of 
short-chain alcohols, methanol and ethanol, on two different fully hydrated lipid bilayer systems in 
the fluid phase at 323 K. It is found that ethanol has a stronger effect on the structural properties 
of the membranes. In particular, the bilayers become more fluid and permeable: Ethanol molecules 
are able to penetrate through the membrane in typical time scales of about 200 ns whereas for 
methanol that time scale is considerably longer, at least of the order of microseconds. We find good 
agreement with NMR and micropipette studies. We have also measured partitioning coefficients 
and the rate of crossing events for alcohols, i. e., typical time scale it takes for a molecule to cross the 
lipid bilayer and to move from one leafiet to the other. For structural properties, two-dimensional 
centre of mass radial-distribution functions indicate the possibility for quasi long-range order for 
ethanol-ethanol correlations in contrast to liquid-like behaviour for all other combinations. 



I. INTRODUCTION 

It is well known that even small changes in the com- 
position of cell membranes can strongly affect the func- 
tioning of intrinsic membrane proteins, such as ion and 
water channels, whi ch regulate the chemical and physi - 
cal balance in cells l|Cantoil l2003t iMazzeo et all Il988(l . 
Such changes may occur due to the introduction of short- 
chain alcohols, or other anaesthetics, at membrane sur- 
faces. Although anaesthetics are being used every single 
day in hospitals around the world, the molecular level 
mec hanisms of general anaesthesia remain elusive, se e 
e.g. (ICantoil[l997tlEckenhofaf20niUTang and Xdl20n2^ . 
The same applies to the effect of alcohols on biological 
systems. iKlemml l)l998|) provides a good review of the 
topic. 

Another aspect to the effect of alcohols appears in a 
more applied context. In the process of producing al- 
coholic beverages, wine in particular, yeasts like saccha- 
romyces cervisiae have to sustain high ethanol concen- 
trations without losing their viability. However, in about 
10% of all wine fermentations the industry encounters 
so-called stuck fermentations (|Bisson and Block! 120021 



Ida Silveira et al.l l2003f) . There is no satisfactory un- 
derstanding of this effect. Some models propose that 
an effect very similar to general anaesthe sia is responsi- 
ble fo r rendering the yeast cells dormant (|Cramer et al.l 
I2002D . It has been suggested that high alcohol concen- 
trations change the membrane structure and force trans- 
membrane proteins into unfavourable conformations. In 
these conformations proteins cannot fulfil their functions 
and thus the yield drops dramatically. 

In addition to the above aspects, there are other 
important issues as well. In particular, in cellular 
systems such as bacteria and yeast, the toxicity of 
ethanol has been suggested to be due to its interaction 
with membranes (|da Silveira et"al l l2fMliT5ri il l2ool 
iLv and Longol l2004(l and the consequent general effects 
such as changes in mechanical properties, permeability 
and diffusion. Changes in such generic membrane proper- 
ties may affect the functions of prot eins and bindingi sites 
due to changes in lateral pressure ijEckenhofl l200lt) . or, 
if the membrane becomes more permeable, changes in 
the electrostatic potential may affect signalling. These 
effects are not to be mixed up with the toxicity due to 
metabolic products such as acetaldehyde from consump- 



2 



tion of ethanol - the cause of poisoning commonly known 
as hangover. 

We concentrate on the effects of ethanol and methanol 
on structural properties of membranes. It is quite surpris- 
ing that despite a vast number of clinical and biochemical 
studies, there has been very few computational investiga- 
tions of the effect of short-chain alcohols, or other anaes- 
thetics, on membranes. The only simulational studies of 
bilayer s and ethanol are, to the authors' knowledge, the 
one bv lFeller et al.l 1^2002^ who used molecular dynamics 
simulations of ethanol and POPC (palmitoyloleoylphos- 
phatidylcholine) lipid bilayers and NMR to study the 
molecul ar level interactio ns in these systems, and the ar- 
ticle by iLee et all ()2004|) discussing alcohol-membrane 
systems briefly. Direct comparison of our results with 
Feller et al. is not meaningful since their study was per- 
formed using a different ensemble close to the gel state 
at 283 K whereas here we are in the biologically relevant 
fluid phase at 323 K. For methanol-bilayer systems there 
exists t o the authors' knowledg e only one computational 
article (|Bemporad et all 1200^ . 

For anaesthetics the situation is slightly better. 
iTanEf and XiJ ((20^3) used molecular dynamics simula- 
tions to study molecular level mechanisms of general 
anaesthesia using halothane as a specific anaesthetic. 
They concluded that the global effects of anaesthetics, 
i.e., due to generic interaction mechanisms, are impor- 
tant and lead to modulations in the functions of channels 
and / or proteins. These conclusions are also supported 
by the fact that the same anaesthetics are effective for 
humans and a variety of animals. Similar conclusions for 
halotha ne interactio n s with bilay ers have been pointed 
out bv iKoubi etldl l|200li |2000() . The importance of 
generic effects has also been indicated in recent exper- 
imental st udies of the e ffect of ethanol on Oenococcus 
oeni cells l)da Silveira et al.^ .2003if) . Although the short- 
age of simulational studies may be due to the high com- 
putational demands of these systems, it is still surprising 
since computer simulations can provide detailed informa- 
tion about fundamental molecular level mechanisms. 

In this article we study the effect of two short-chain al- 
cohols, ethanol and methanol, on two different lipid mem- 
branes consisting of either pure DPPC (dipalmitoylphos- 
phatidylcholine) or POPC. Methanol is a small so- 
lute having a single hydrophilic hydroxyl group whereas 
ethanol possesses an additional hydrophobic carboxyl 
group. DPPC and POPC share the same headgroup 
but one of the tails of POPC has a double bond and is 
two carbon atoms longer, whereas DPPC has only single 
bonds in its chains, see Fig.Q] We have studied these sys- 
tems under fully hydrated conditions using microscopic 
molecular dynamics. 50 ns trajectories for each of the 
four combinations of lipid and alcohol allow us to gather 
high statistical accuracy. 

Phospholipid bilayers can be considered as a first ap- 
proximation to understand the behaviour of cell mem- 
branes under the influence of alcohol, and much infor- 
mation can be extracted from such systems. The simula- 



tions show that ethanol is able to pass through the bilayer 
much more easily than methanol. This can be explained 
by the hydrophobic nature of the carbon "tail" of ethanol, 
making passing through the hydrophobic tail regions 
of lipid bilayers easier. In addition, ethanol molecules 
condense near the interface region between lipids and 
the surrounding water, i.e., there is a sharply increased 
density of ethanol near the interface region, while for 
methanol only a moderate increase of the density is seen 
near the interface region. These effects are very pro- 
nounced for DPPC bilayers, and only slightly weaker for 
POPC bilayers. This has far reaching implications for 
the basic properties of bilayers. 

The rest of this article is organised as follows. In the 
next section we describe the model and the simulation 
details. Then, in Sec. IIIII we present the results from 
the simulations. Section IIVI contains a discussion and 
conclusions. 



II. MODEL AND SIMULATION DETAILS 

We have simulated lipid bilayer systems consisting of ei- 
ther 128 DPPC or 128 POPC molecules (i.e., 64 lipids in 
each leaflet). For th e lipids we used a previously va lidated 
united atom model (iTieleman and Berendsenl 1199611 . The 
DPPC simulations are based on the final structure of a 100 ns 
run of a DPPC bilayer that is fully hydrated by 3655 water 
molecules . The configuratio n is available online^. The simu- 
lations bv lPatra eTayillooi) were run using the same param- 
eters as here (details below), and the 100 ns run used i n this 
study is a continu ation of a 50 ns study iFalck et al.L 120041 : 
iPatra et al^.l2004^ . For the POPC simulations, such an initial 
structure had to be gen erated first. We start ed with a fully 
hydrated POPC bilayer jTieleman et al.Ul999ll and simulated 
it for 10 ns. The final structure of that simulation run was 
used as a starting point for the POPC simulations reported 
here. 

In order to add the alcohol molecules, the simulation box 
was first extended in z-direction such that an empty vol- 
ume was created. A total of 90 ethanol (methanol) molecules 
were randomly inserted in the empty volume, and the remain- 
ing space was filled with water. The total number of water 
molecules then amounted to 8958 for the DPPC systems and 
8948 for the POPC systems, or, in other words, 1 mol% alco- 
hol on a lipid free basis. The small difference between DPPC 
and POPC systems is due to the different lateral extensions 
of the bilayers. 

The force field parameters for bonde d and non- 
bonde d interaction were taken from iBerger et alJ 
il997f). available online^. Partial charges were taken 
from ITieleman and Berendsenl il996t) . available online for 
both DPPC ^ and POPC As is seen in the chemical 
structures in Fig. [H DPPC and POPC are identical up to 
a single pair of CH-groups, connected by a double bond 



^ http:/ /www. softsimu.org/downloads.shtml 
^ http:/ /moose. bio. ucalgary.ca/Downloads/filcs/lipid.itp 
^ http:/ /moose. bio. ucalgary.ca/Downloads/files/dppc.itp 
http:/ /moose. bio. ucalgary.ca/Downloads/files/popc.itp 
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10 water molecules per lipid. This was done in order to see 
if dehydration has a direct effect on the properties but no 
significant effects were found. 

The simulations took a total of 16 000 CPU hours using an 
IBM eServer Cluster 1600 (Power4 processors). 



III. SIMULATION RESULTS 
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FIG. 1 Structures of POPC (top) and DPPC (bottom). They 
are identical with the exception of the sn-2 chain which is 
two carbons longer and contains one double bond for POPC. 



instead of a single bond in the sn-2 chain of POPC, and the 
two additional CH2 groups at the end of that chain. This 
similarity is reflected in the force flelds, which are identical 
up to the modelling of the four affected atoms. Ethanol 
and methan ol were modell e d usi ng the Gromacs force field 
parameters (iLindahl et alj. I2OOIII which are identical with 
the exception of the added CH2 group for ethanol. Thus, 
differences observed between the two lipids or the two alco- 
hols do not originate from differences in their respective force 
field parameterisations but are due to the physics and / or 
chemistry of those co mponents. For water th e Simple Point 
Charge (SPC) model iBerendsen et allll98lfl was used. 

T he simulations were performed with the Gromacs pack- 
age iLindahl et al.ll200lfl . The lipids, water molecules and al- 
cohols were separately coupled to a heat bath at tem perature 
T = 323 K using the Berendsen thermostat (iBercn d scn et alJ . 
Il984fl with a coupling time constant of 0.1 ps. All the bond 
lengths were con strained to their equilibrium values by the 
Lines algorithm (iHess et al.L 1199711. Pressure was co ntrolled 
using the Berendsen barostat iBerendsen et all . 1198411 with a 
time constant of 1 ps. The pressure coupling was used semi- 
isotropically such that height of the box {z direction) and 
the cross sectional area (xy-plane) were allowed to vary in- 
dependently of each other. Lennard- Jones interactions were 
cut off at a distance of 1.0 nm and the time step was set to 
2 fs. Long-range electrostatic s were updated every 10-th time 
step [the twin-range scheme iKessel et alJ . l200^ iPatra et all 
l2003f) was used], and handled b y the particle-mesh Ewald 
(PME) algorithm iEssman et all Il995ll . For DPPC bilay- 
ers it has been shown that replacing PME by the compu- 
tation ally cheaper cut-off sch eme leads to pronounced arti- 
facts (iPatra et all l200l I200I) . 

The systems were simulated for a total of 50 ns each. After 
20 ns, the samples had equilibrated, and the remaining 30 ns 
were used for data collection. Equilibration was determined 
by monitoring the area per lipid as described in the next sec- 
tion. For completeness, we also present results for pure DPPC 
and POPC bilayers. In particular the latter ones are impor- 
tant since many of the results have so far only been published 
based on simulations using a cutoff for handling electrostat- 
ics. In addition to the above systems, we also performed a 
control simulation with a dehydrated system containing only 



Before presenting a systematic summary of our results, 
we give a quick overview of the basic properties of these 
systems. Alcohol molecules have a tendency to collect in 
or near the bilaycr (Scc. lIII.B|) . This tendency is stronger 
for ethanol than for methanol as confirmed by a partition 
analysis (Sec. IIII.III . Ethanol is able to form hydrogen 
bonds with the lipids in the bilayer fScc. lIII.D|l . and these 
hydrogen bonds reduce the order parameter of the lipid 
hydrocarbon tails. The combination of all this results 
in an easy penetration of ethanol through the bilayer. In 
contrast, no hydrogen bonds or penetration was observed 
for methanol. 

In this paper we use the following colour code for all 
figures. Curves for systems containing ethanol arc drawn 
in red, curves for methanol in green and pure lipid sys- 
tems without alcohol in blue. 



A. System dimensions 

The area per lipid is one of the most important quan- 
tities characterising lipid bilayer systems and it can also 
be used to monitor equilibration during a simulation run. 
The time evolutions of the area per lipid in the systems 
studied here are shown in Fig. [3 The average areas per 
lipid, (A), obtained in our simulations are listed in Ta- 
ble |I| 

For a pure DPPC bilayer we obtain (Adppc) = 
0.655 nm^ agreeing well with previo u s sim ulations and 
experiments, see Ref. (|Patra et all 1200,*^ and refer- 
ences therein. For pure POPC we obtain (^popc) = 
0.677 nm^ in agreement with previous computational 




FIG. 2 Temporal behaviour of the area per lipid A{t) for a 
DPPC bilayer (left) and a POPC bflayer (right). The colour 
of the line marks whether the lipid bilayer has been simulated 
in the presence of ethanol (red), of methanol (green), or of no 
alcohol (blue). This is our standard colour code employed 
throughout this paper. 
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System 


Average area per lipid 


DPPC (pure) 


(0.655 


± 


0.002) 


nm^ 


DPPC + ethanol 


(0.699 


± 


0.002) 


nm^ 


DPPC + methanol 


(0.693 


± 


0.004) 


nm^ 


POPC (pure) 


(0.677 


± 


0.003) 


nm^ 


POPC + ethanol 


(0.699 


± 


0.003) 


nm^ 


POPC + methanol 


(0.693 


± 


0.003) 


nm^ 



TABLE I Average area per lipid for all systems studied in 
this work. A weak effect of the alcohols is visible. The er- 
ror estimates have been computed from block averaging and 
extrapolating to large block sizes. 



studi es ()Chiu et all Il999t iPasenkiewicz-Gierula et all 
|2003|) and slightly larger than the resu lts from x-ray 
diffraction studies l|Pabst et all l2000alhD . For POPC, 
the difference to x-ray diffraction results (Apopc) ~ 
0.61 nm^ may be due to differences in trans-gauche con- 
formational changes. 

As seen from Table the presence of alcohol has a 
small but non-vanishing effect on the area per lipid. The 
number of water molecules per lipid molecule plays only 
a minor role as was verified by an additional simulation 
of DPPC with ethanol and a reduced amount of water. 
Interestingly, ethanol and methanol have almost the same 
effect on the area per lipid. 

While the definition of the area per lipid is straight- 
forward, the same is not true for the volume occupied by 
a lipid. The precise definition of the volume V (or the 
thickness d) of a memb r ane is non-trivial as discussed at 
length bv lArmen et alJ l)l998() . Here, we chose an opera- 
tional definition based on local mass density. Other defi- 
nitions, e.g., employing the electron density, are equally 
possible. 

Below, we give the two definitions we used to compute 
the thickness. If piipid, Pwatcr and Paicohoi are the mass 
densities of the three components, the effective thickness 



System 


di [nm] 


Vi [nm^] 


d2 [nm] 


V2 [nm^] 


DPPC (pure) 
with ethanol 
with methanol 


2.02 ±0.05 
1.84 ±0.02 
1.93 ±0.06 


1.32 ±0.03 
1.28 ±0.02 
1.34 ±0.04 


2.02 ±0.05 
1.90 ±0.02 
1.95 ±0.06 


1.32 ±0.03 

1.33 ±0.02 
1.35 ±0.04 


POPC (pure) 
with ethanol 
with methanol 


1.96 ±0.04 
1.87 ±0.02 
1.94 ±0.02 


1.33 ±0.03 
1.31 ±0.02 
1.35 ±0.01 


1.96 ±0.04 
1.93 ±0.03 
1.96 ±0.02 


1.33 ±0.03 

1.35 ±0.02 

1.36 ±0.01 



TABLE II The thickness d of a leaflet (the bilayer thickness 
is twice that value) and the corresponding volume per lipid 
using the two definitions given in Eq. Q. The error estimate 
for d has been computed by cutting the analysis part of the 
trajectory in two parts and by applying Eq. Q separately to 
both of them. 



of a single leaflet can be defined by 

di = - dz , (la) 

2 J Plipid + Pwatcr + Palcohol 

d2 = U '-^ dz . (lb) 

^ J Plipid ' Pwatcr 

These two definitions differ in their treatment of the 
alcohol volume fraction and give the same thickness for 
pure lipid bilayers. After defining the thickness, the vol- 
ume is simply V = d{A) with (A) being the average area 
per lipid. The results using both of the above definitions 
are summarised in Table ITll 

The thicknesses we obtained for pure POPC agree 
very well with recent x-ray diffraction studies of 
Vogel et al.| (| 200Cl|) and c o mpute r simulation studies of 
Gullingsrud and SchulteiJ l)2004(l who obtained 3.9 nm 
and 3.92 nm, respectively, for the total bilayer thick- 
ness 2 d. For DPPC the thickness and volume 
are a few percent larger than the experimental re- 
sults (|Nagle and Tristram-Nad j. l2000|) . Using the elec- 
tron density to define the thickness would have led to 
similar results, see Sec. IIII.Cl 

A comparison of Tables and |n] shows that the ad- 
dition of ethanol or methanol to a bilayer expands its 
surface slightly while the thickness decreases such that 
the volume per lipid does not change significantly. This 
is as assumed since the main effect of the addition of al- 
cohol is a reduction of the surface tension of the water 
phase. This is in agreem ent with observ ations from a 
DPPC-halothane system (iTu et alJ.ll998|) . We will re- 
turn to this issue in Sec. II VI 

To complement the average area per lipid (A) mea- 
surements, we have also computed the area probability 
distribution P{A) by Voronoi tessellation. By definition, 
Voronoi tessellation measures the area that is closer to 
a given molecule than to any other one. The Voronoi 
approach does not uniquely specify which point should 
be used to represent the entire molecule. We used the 
centre-of-mass position of each lipid, projected onto the 
xy-plane. Other choices are also possible, such as the po- 
sition of the sn - 3 carbon which gives a better indication 
of the backbone of the lipid whereas the centre-of-mass 
describes the entire lipid. 
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FIG. 3 Distribution of the individual areas of the lipids as de- 
termined by two-dimensional Voronoi tessellation for DPPC 
(left) and POPC (right). 
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FIG. 4 Autocorrelation function for the individual areas of 
the lipids as determined by two-dimensional Voronoi tessella- 
tion for DPPC (left) and POPC (right). 
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FIG. 6 Mass density profiles across the bilayer for DPPC 
(left) and POPC (right). The density of the alcohol has been 
scaled by a factor of 10, and the colour code is the same as 
above. 



The resulting distributions are shown in Fig. 

The mean of that distribution is, by construction, identi- 
cal to the average area per lipid as shown in Tabled and 
thus does not contain any additional information. For 
this reason, not the plain distribution P(^A) but rather 
the re-scaled distribution is shown in Fig. O 

Plotting the result in this way shows that alcohol does 
not influence the Voronoi distribution in any way that 
would not be captured already by the average area per 
lipid. 

It is also possible to compute the autocorrelation time 
of the Voronoi areas. This time gives an indication of 
how quickly the geometry of the bilayer changes locally 
whereas the fluctuations in the size of the simulation box 
seen in Fig. |2] are related to global changes of the geome- 
try. The results are shown in Fig. 21 The faster decay for 
the systems with ethanol suggests that the bilayer might 
become more fluid but care should be taken in draw- 
ing conclusions as the differences between the curves are 
rather small. 

Next, we perform similar Voronoi tessellation for the 
alcohol molecules inside the bilayer interface region. The 
precise definition of that region turned out not to be crit- 
ical, and we included all alcohol molecules within the 
range 0.7 nm < z < 1.7 nm from the centre of the bi- 
layer. (Our choice for this range is motivated by the 
results to be discussed later in Sec. IIII.Bl ) The variable 
number of molecules forbids a proper calculation of the 



correlation time for the areas assigned to each alcohol, 
though, and thus wc present only the distribution P{A) 
in Fig. [SI Since there are fewer methanol molecules close 
to the bilayer than there are ethanols (sec Fig.O, the av- 
erage area per methanol is larger than the average area 
per ethanol. 



B. Mass density 

The mass density profiles across the bilayer are shown 
in Fig. 1^1 For each analysed simulation frame, the system 
was first translated such that the centre of the bilayer 
was located at z = 0. Particles with z < were mirrored 
to z > to reduce statistical error. The masses of the 
hydrogen atoms were accounted for in the calculation. 
Due to the low density of alcohol, its curve is scaled by 
a factor of 10 in the figure. 

Additional information can be gained by considering 
separately the two charged groups in the lipid head- 
groups, namely the phosphate (P) and the choline (N) 
group (cf. Fig.P). In addition to this, the oxygen atom of 
the alcohols is included in Fig.[7| We could not find a di- 
rect comparison for the mass densi ty but the observat ions 
from the computer simulations of lFeller et all 1I2OO2D are 
consistent with our results for the mass density. Fig- 
ure [7| shows that the alcohol molecules have a strong 
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FIG. 5 Distribution of the Voronoi areas of the alcohols with 
DPPC (left:) and POPC (right). Only alcohol molecules lo- 
cated close to the bilayer interface are included in the analysis. 
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FIG. 7 Mass density profiles across the bilayer. The densities 
given are the O(xygen) of alcohol as well as N(itrogen) and 
P(hosphorus) of the lipid, and scaled to give a maximum of 
unity. Left for DPPC, right for POPC. 
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FIG. 8 Density profile across the whole bilayer. The lipid 
component is divided into the contributions from the two sep- 
arate leaflets, providing a measure of interdigitation. The al- 
cohol component has been suppressed in the figure. Left for 
DPPC, right for POPC. 

tendency to accumulate below the bilayer-water inter- 
face layer (approximately given by the location of the 
phosphate and choline groups) , and that this tendency is 
stronger for ethanol than for methanol. We will return 
to this issue in the partitioning analysis in Sec. IIILll and 
the membrane penetration analysis in Sec. IIII.JI 

The density of the lipid is decreased near the centre of 
the bilayer. This phenomenon is known as lipid trough 
and means that the two leaflets are repelling each other. 
Still, the tails of the lipids from one leaflet are able to 
penetrate into the other leaflet; thi s is known as inter- 
digitation (jLobbcke and Cev3 . ll995t) . To analyse this, in 
Fig. IHlwe plot the density throughout the whole bilayer, 
i. e., the positions of all atoms are not folded into a single 
leaflet. We show separately the density of the lipids be- 
longing to the upper and the lower leaflet of the bilayer. 
It is easily seen that the tails of the lipids can penetrate 
up to 0.5 nm into the other leaflet, and the degree of 
interdigitation is largely independent of the presence of 
alcohol. 



C. Electron density 

Electron densities provide information about the struc- 
ture of bilayers along the normal to the bilayer plane sim- 
ilar to the mass densities. Experimentally, x-ray diffrac- 
tion provides a means to access this quantity, the mea- 
surements yielding information of the total electron den- 
sity profile. 

Figure shows the total electron densities in different 
cases. For the pure lipid bilayers, the curves show the 
typical behaviour, namely a maximum associated with 
the electron dense areas in the headgroup, i. e., the phos- 
phate groups, and the mini mum at the bilayer centre - 
the s o-called methyl trough l)Nagle and Tristram-NaglS 
|2000() . Experimentally, an electron density profile con- 
tains different information than a mass density profile 
since the chemical composition at depth z is not known 
directly. (For computer simulations this problem does 
not exist.) 

For the pure systems our results agree well with ex- 



FIG. 9 Electron density profiles in the studied systems with 
DPPC (left) and POPC (right). 

periments (|Nagle et al.Lll99^ . The only x-ray diffraction 
study in a related system conta ining ethanol o r methanol 
that we are aware of is that of lAdachil l)2000l) . Unfortu- 
nately a direct comparison is not meaningful since that 
study was done in the gel phase using multilamellar vesi- 
cles as compared to the fluid-like phase and planar bilayer 
system studied here. 



D. Hydrogen bonding of alcohol to lipids 

As was shown in Sec. lIII.Bl the alcohol molecules have 
a tendency to be located inside the bilayer, and this ten- 
dency is stronger for ethanol than for methanol. The 
alcohol molecules are not located directly at the water- 
membrane interface but rather further inside the bilayer. 
For the simulations with ethanol, a direct visual inspec- 
tion of the atom positions shows that ethanol molecules 
are located close to the ester oxygens of the lipids, see 
Fig.Uni 

This visual conclusion is confirmed by a hydrogen 
bonding analysis. In such an analysis, possible donors 
and acceptors are identified by their chemical properties, 
and a hydrogen bond is then assumed to exist whenever 
two such atoms and an additional hydrogen atom ful- 
fil certain geometric conditions. (The distance between a 
hydrogen and an acceptor has to be smaller than 0.25 nm, 
and the angle between acceptor, hydrogen and donor has 
to be smaller than 60 degrees.) 

The hydrogen bonding analysis shows that the major- 
ity of the ethanols are involved in hydrogen bonds with 
lipids whereas not a single hydrogen bond between a 





DPPC POPC 


bound ethanols 


72.9 


71.6 


bound lipids 


59.7 


59.2 


hydrogen bonds 


74.1 


72.8 


lifetime [ns] 


1.20 


1.15 



TABLE III Results of the hydrogen bonding analysis for the 
DPPC and POPC bilayers with ethanol. (The systems consist 
of 128 lipid molecules and 90 alcohol molecules.) No hydrogen 
bonds between methanol and lipids were found. 
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FIG. 10 A DPPC molecule together with two ethanol 
molecules. The ethanols are located close to the ester oxygen. 
The DPPC molecule is drawn as rods whereas the ethanols 
are drawn in a spacefilling representation. To aid the eye, the 
ethanols are coloured blue-ish. 



methanol and a lipid molecule was found in our simu- 
lations. The results are summarised in Table Hill Many 
lipids are involved in more than one hydrogen bond which 
comes as no surprise since they possess an ester oxy- 
gen in each of their two chains. Comparison with the 
lifetim e data for ethanol in Table IIIII with NMR experi- 
ments (|Holte and Gawrischlll99^ shows excellent agree- 
ment. In their experiments Holte and Gawrisch reported 
the lifetimes to be around 1 ns while we obtained 1.20 ns 
for the ethanol lipid hydrogen bonds. We are not aware 
of any such experiments for methanol. 

The number of alcohol molecules involved in hydro- 
gen bonds is best compared against the total number of 
alcohol molecules located inside the bilayer. The latter 
number is relatively ill-defined but from Fig. one can 
easily compute that for cthanol-containing systems only 
of the order of 10 ethanol molecules out of the approxi- 
mately 70 inside the bilayer are not involved in hydrogen 
bonds. For comparison, in the methanol systems there 
are of the order of 20 methanol molecules inside the bi- 
layer and none of them is involved in hydrogen bonds (but 
cf. Sec. IIII.Jll . These numbers show that there indeed is 
a significant difference between ethanol and methanol. 

Hydrogen bonding analysis offers a well-defined crite- 
rion to decide whether a given lipid is interacting strongly 
with an alcohol molecule or not. This will be used in the 
following sections to study separately the two lipid pop- 
ulations, lipids bound to an alcohol and lipids not bound 
to an alcohol. 



E. Radial-distribution functions 

In addition to the mass density profiles, valuable infor- 
mation may be gained from radial-distribution functions 
(RDFs). The radial-distribution functions g{r) give the 
probability of finding two particles at a mutual distance r 
once geometric and density factors have been scaled out. 

Figure ^2 shows the RDFs between the oxygen of the 
alcohol and different charged groups inside the lipid. 
Some of these groups were depicted already in the mass 
density profile in Fig. but whereas there only the verti- 
cal distance between particles was considered, the RDF 
considers the real three-dimensional distance between 
them. 

While the mass density profile showed that on the av- 
erage ethanol molecules prefer to reside 0.5 nm below 
the lipid headgroups, the RDF shows that the three- 
dimensional preferred distance is only 0.38 nm. This is 
no contradiction but is easily understood by the observa- 
tion (cf . Fig. El £^ bit further down) that lipid molecules 
without an attached ethanol molecule are sticking out of 
the bilayer more than those with an attached ethanol. 
This is captured only by the RDFs but not by the mass 
density profile. 

The radial-distribution functions for the different sys- 
tems look quite similar - with one exception: The RDF 
between the alcohol and the ester group is peaked at a 
much smaller distance for ethanol than it is for methanol. 
This is in agreement with the results of the hydrogen 
bonding analysis. 

By studying the mutual RDFs of the choline and / or 
phosphate groups, it is possible to detect phase transi- 
tions of the bilayer. Within error margins, these RDFs 
are not dependent on the presence of alcohol, and for 
space reasons we do not show them here as they are 
identical to the RDFs published previously ijPatra et al.l 
I2004D . 

We have also studied two-dimensional radial- 
distribution functions of entire molecules, i.e., the 
molecules' centres-of-masses were projected onto the 



N(lipid)-0(alco.) 
P(lipid)-0(alco.) 
O(lipid)-0(alco.) 




-i/4 


I N(lipid)-0{alco.) - 




1 P(lipid)-0(alco.) --- _ 




\ 0(lipid)-0(alco.) 


J : 1 








. 1 . i" : •■i 



2 0.4 0.6 0.8 1.0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



r [nmj 



r [nmJ 



FIG. 11 Radial-distribution function between the oxygen of 
alcohol, on the one hand, and the phosphorus and the nitrogen 
atoms in the headgroup as well as the ester oxygens in the 
lipid tail, on the other hand. The oxygen-oxygen curve has 
been scaled by a factor of 1/4, i.e., in reality the RDF peaks 
at a value four times as large as displayed in the figure. 
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FIG. 12 Two-dimensional radial-distribution functions of the centre-of-mass positions for lipid-lipid (left), lipid-alcohol (centre) 
and alcohol-alcohol (right). 



xj/-plane and radial-distribution functions were then 
computed. The results are shown in Fig. 1121 The mutual 
radial-distribution functions of the lipids exhibit a very 
soft core as lipids arc able to wrap around each other. 
No dependence on lipid type or presence of alcohol was 
observed. 

The RDF between alcohol and lipid is qualitatively dif- 
ferent for ethanol and methanol. For an ethanol, there is 
a large probability for it to be at the same x-y position 
as the (centre-of-mass of the) lipid. This reflects the hy- 
drogen bonding of ethanol close to the centre of the lipid. 
This bonding is absent for methanol, and consequently 
then g{r) for r — > 0. No significant dependence on 
the kind of lipid is observed. 

The alcohol-alcohol radial-distribution functions have 
a very different character: For methanol with DPPC, 
the first peak is very distinct but the correlations de- 
cay soon after. For methanol with POPC an additional 
peak is observed. (This is the only curve with a signif- 
icant difference between DPPC and POPC. We cannot 
offer a convincing explanation for this.) For ethanol, the 
behaviour shows almost quasi long-range order. The rea- 
son for this ordering is not clear and further experiments 
would be needed to study this in detail. 



F. Order parameters 

Ordering of the lipid acyl chains is typically charac- 
terised using the order parameter tensor 

Sap = ^ {ScOsOa COS 6,3 - S^p) , (2) 

where a, (3 = x, y, z, and 9a is the angle between the a'^ 
molecular axis and the bilayer normal (z-axis). The order 
parameter is then computed separately for all carbons 
along the acyl chain. Since lipid bilayer systems possess 
symmetry with respect to rotations around the z-axis, 
the relevant order parameter is the diagonal clement Szz ■ 
To relate Szz to the experimentally relevant deuterium 
order parameter 

2 1 

•S'cD = -^Sxx + -^Syy , (3) 

we use the symmetry and write Sxx = Syy^ and Sxx + 
Syy + Szz = 0. Using these relations we have S'cD = 
—Szz/'^- To allow comparison with experimental data, 
we present our results in terms of I^cdI- 

Since our simulations employ a united atom model, 
no explicit information about the hydrogen positions is 
available and they must be reconstructed assuming a per- 
fect tetrahcdral arrangement. The results are shown in 
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carbon atom 

FIG. 13 Deuterium order parameters, computed from Eq. I^, 
for DPPC (solid line) and POPC (dashed line). For DPPC, 
the average over the two tail chains is displayed while for 
POPC only the saturated sn-1 chain is shown. For number- 
ing of carbon atoms, refer to Fig. Q 
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FIG. 14 Order parameter of DPPC (left) and POPC (right) 
in the presence of ethanol. For DPPC, both the sn-1 and 
sn-2 chains are shown whereas for POPC only the sn-1 
chain is depicted. The order parameter has been computed 
separately for lipids that are bound to at least one alcohol 
molecule (solid lines), and for lipids that are not bound to 
any alcohol (dashed lines). 
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FIG. 15 Distribution of the angle 6 between the P-N vector 
and the bilayer normal for DPPC (left) and POPC (right). 

Fig. El Reconstruction of the hydrogen positions means 
that for the outermost carbon atoms of the tail no order 
parameter can be constructed. This also includes po- 
sitions where a sequence of carbon atoms connected by 
single bonds ends in a carbon atom having a double bond. 
It is not a problem to compute the order parameter for 
a chain of atoms connected by double bonds, there just 
is a problem connecting such a chain to a chain of atoms 
connected by single bonds. For this reason we show the 
order parameter for both chains of DPPC but for POPC 
we restrict ourselves to the saturated sn-\ chain, sec 

Fig.m 

The results for the order parameter for all the cases are 
shown in Fig. El We find that methanol increases order- 
ing of the lipid acyl chains close to the glycerol group, 
while the effect of ethanol is strongest below the glycerol 
group, around the centre of the hydrocarbon tails. These 
results are fully consistent with the mass density profiles 
in Fig. El 

As discussed in the introduction, there is, to oxvc knowl- 
edge, only one other, related membrane-alcohol simula- 
tion study l)Feller et al.l l2002() . In that the temperature 
was 40 degrees lower. Although direct comparison is not 
possible, the order parameters in that study and ours are 
in qualitative agreement. It is noticeable that the effect 
of alcohols on the average order parameter is small - this 
is consistent with the fact that the volume occupied by 
a lipid changes only slightly. 
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FIG. 16 Distribution of the angle between the P-N vector 
and the bilayer normal for DPPC (left) and POPC (right) in 
the presence of ethanol. The distribution is separated into 
lipids that are hydrogen bonded respectively not hydrogen 
bonded to an ethanol molecule. 




FIG. 17 Distribution of the position of the phosphate group 
(P), the choline group (N) and the centre-of-mass of the en- 
tire lipid (COM), left for DPPC, right for POPC. Results are 
shown as solid (dashed) lines for lipid molecules bound (un- 
bound) to an ethanol molecule. Please note that in this figure 
the colours do not mark the kind of alcohol present. 



Additional insight can be gained by combining the 
order parameters with the hydrogen bonding analysis 
from Sec. IIII.DI This combination allows the compu- 
tation of the order parameter depending on whether the 
lipid forms a hydrogen bond with an alcohol molecule or 
not. The result in Fig. 1141 shows that binding of ethanol 
slightly decreases the order of the tails. 

In addition to the deuterium order parameter for the 
acyl chains, it is possible to study the ordering of head- 
groups in a similar fashion. To do that, we have chosen 
the angle of the P-N vector with respect to the bilayer 
interface plane and have computed its distribution. The 
result is shown in Fig. El The P-N vector has a sig- 
nificantly higher tendency of being in the bilayer plane 
{9 = 90°) than of sticking out of it. The computed dis- 
tribution is only slightly dependent on the presence of 
alcohol. 

Again, we can obtain additional information for the 
ethanol systems if the angular distribution is separated 
into the distributions of lipids that are hydrogen bonded 
to an alcohol, and those that are not. The result is shown 
in Fig. El For DPPC the angular distribution is not 
infiuenced at all by hydrogen bonding whereas for POPC 
the influence is small. 

To better explain the results presented in this section, 
we return to the mass density profiles. We analysed the 
positions of the phosphate and the choline groups in the 
heads of the lipids as well as the lipids' centre-of-mass po- 
sitions, separated into lipids that are hydrogen bonded or 
not hydrogen bonded to an ethanol. The data in Fig. 1171 
shows that lipid molecules are shifted towards the centre 
of the bilayer by approximately 0.2 nm if they are bonded 
to an ethanol molecule. This is in agreement with the 
reduced order parameter as that is normally associated 
with a thinner bilayer. The vertical distance between the 
choline and the phosphate group remains unchanged, in 
agreement with the distribution of the P-N angle. 



10 



0.05 

0.00 
-0.05 
-0.10 
-0.15 
-0.20 
-0.25 
-0.30 
-0.35 
-0.40 



.1 1 1 1 1 1 


1 i i 1 1 1 1 1 1 1 i i 1 1 1 1 L 










1 . . . . 


\y POPC - i 

.... 1 .... 1 .... 1 . ; 



1.0 1.5 



2.0 

Z [nm] 



2.5 3.0 



FIG. 18 Water orientation, as described by the mean cosine 
of the angle of the water dipole moment with respect to the 
bilayer normal. 



G. Orientation of the water dipole 

Ordering of the water dipole in the vicinity of the 
bilayer-water interface is described by calculating the 
time averaged projection of the water dipole unit vector 
fl{z) onto the interfacial normal n, 



P{z) {fl{z) ■ n) = (cos6l(z)) , 



(4) 



where z is the 2-component of the centre-of-mass of the 
water molecule and vector n points away from the bilayer 
centre along the ^-coordinate. 

The data for all the studied cases are shown in Fig. 1181 
For pure bilayers the results are i n agre ement with pre- 
vious studies, e.g. (|Patra et al.L l2003j) . When either 
methanol or ethanol is added, the water dipole becomes 
less oriented, i. e., the addition of alcohol slightly reduces 
the amount of ordering. For pure bilayers, and for bilay- 
ers with added methanol, the minimum remains at the 
same distance, at about 1.8 nm from the bilayer centre. 
For added ethanol the minimum shifts to a smaller dis- 
tance, to about 1.6 nm. This is a reflection of the fact 
that ethanol leads to a slightly larger area per lipid and 
thus a thinner bilayer (see Sec. IIII.^ . 



H. Electrostatic potential 

To obtain the electrostatic potential across the bilayer 
the average charge density profile was first computed such 
that the centre of the bilayer was separately aligned to 
2 = for each simulation frame. Then, the electrostatic 
potential was determined by integrating the charge den- 
sity twice with the initial condition V{z = 0) = 0. 

The electrostatic potentials for all studied cases are 
shown in Fig. ^1 For pure DPPC the electro- 
static potential was determined t o be —589 mV in 
agreement with previous studies ijPatra et al.l l2003t 
Iriclcman and Bcrcndscn, 1996). For pure POPC we ob- 
tain -507 mV. 

The addition of alcohol leads only to small changes 
in the electrostatic potential. This is what is expected 
from the results presented so far. The only way in which 
the membrane potential could be changed significantly 




FIG. 19 Electrostatic potential through the bilayer for DPPC 
(top) and POPC (bottom). 



would be by re-arrangement of the P-N angle of the 
headgroup, resulting in a change of the dipole moment 
of the lipid headgroup. No such re-arrangement was ob- 
served, cf. Fig. [T51 

On a superficial level this may seem to be in contra- 
diction to some previous suggestions that the narcotic 
effects of alcohols are mainly due to a change of the elec- 
trostatic potential. We would like to point out, however, 
that even though the direct effect of alcohol to the po- 
tential is small, this does not exclude secondary effects 
which may lead to a significant change in the electrostatic 
potential. We return to this issue in Sec. IIVI 



I. Partitioning coefficients 

To measure whether alcohols prefer to stay close to 
the membrane or rather to be in the water phase, we 
computed the partitioning coefficients. They provide a 
thermodynamic quantity characterising a molecule's ten- 
dency to choose its environment. 



The partitioning coefficient Kp is defined as 

-^'^p — ^bilayer /"^ water : 



(5) 



where .X'hiiavpr and Xwatpr are th e moles of solute per 
kg of solvent l|Westh et all l200l|) . Equation (jsl can be 
reformulated in terms of the molar concentrations n and 
the molar masses m as 



n _ „ water N bilayer 

"-alcohol i "alcohol '"water 



\^ "alcohol/ "alcohol '"lipid 



(6) 



Here, means the molar concentration of component b 
in the phase a. 

One might be tempted to count the number of alcohol 
molecules in the lipid phase and in the water phase us- 
ing some functional definition of where the lipid bilayer 
interface is located. This approach, however, gives incor- 
rect results as the relevant quantity is not the number 
of alcohol molecules inside the membrane but rather the 
number of alcohol molecules influenced by the membrane. 



11 



system 


water 
"'cthanol 


„iipid 

cthanol 




DPPC-ethanol 
POPC-ethanol 


(3.98 ±0.06) • 10"" 
(3.34 ± 0.04) • 10""' 


0.403 ± 0.001 
0.405 ± 0.001 


41.6 ±0.6 
48.3 ±0.9 


DPPC-methanol 
POPC-methanol 


(4.53 ±0.01) • 10^=* 
(4.72 ± 0.02) • 10~3 


0.278 ±0.001 
0.272 ±0.001 


2.07 ±0.02 
1.87 ±0.02 



TABLE IV Results of the partitioning analysis. 



Starting point arc the mass densities Pwatcr and Paicohoi 
far from the bilayer in the water phase. These are avail- 
able from Fig.|51 and knowledge of the molar masses gives 
'T-watcr- Multiplying riwatcr by the total number of water 
molecules in the simulation box gives the number of al- 
cohol molecules that would be there if there was no lipid 
bilayer. The remaining alcohol molecules thus must be 
"brought in" by the lipid bilayer. 

Our results are summarised in Table IIVI Although 
experimental measur ements of partitioning coeflicients 
have been performed llRowe et allll998t iTrandum et all 

WD, 



I2OOOI: IWesth and Tranduml. 11*999(1 . we were not able 
to find experimental data to make exact quantita- 
tive co mparisons. However, the experimental data for 
DMPC l|Trandum et al.U200(]D with cthanol is in quahta- 
tive agreement with our results. They are also in agree- 
ment with the general trends that Kp increases as the 
lipid tails or the alcohols become longer. As we have 
studied the system at single temperature we cannot com- 
ment on the temperature dependence of the partitioning 
coefficients. 



J. Crossing events 

Next, we analyse the penetration of alcohol through 
the membrane. A quick overview can be obtained by 
plotting the z-component of the positions of the alcohol 
molecules. (The positions of all atoms have to be trans- 
lated for every simulation frame such that the centre of 
the bilayer does not move). The result is shown in Fig. 1201 
The density (in space and time) of alcohol molecules is 
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FIG. 20 z-positions of all alcohol molecules as a function of 
time for DPPC (left) and POPC (right). Ethanol (red) is able 
to penetrate into the bilayer (located at z = 0) much better 
than methanol (green). Crossing events of ethanol are seen 
while they are completely absent for methanol. 



large in most parts of the diagram such that it is difficult 
to visually identify individual trajectories. 

It is easily seen from Fig. 1501 that the density of alco- 
hol molecules is reduced in the centre of the bilayer. In 
addition, the gap in the centre of the bilayer is smaller 
for ethanol than it is for methanol. This is in agreement 
with the mass density profiles presented in Sec. IIII.BI It 
is also evident that there is a significant number of events 
where an ethanol molecule is crossing from one leaflet 
to the other while no such events are seen for methanol 
molecules. (Crossing events cannot be inferred from the 
mass density since these events happen so fast that the 
resulting mass density of alcohol in the centre of the bi- 
layer is negligible.) 

In the simulations, it is directly known which atoms 
form the lipid molecules of the upper leaflet of the bi- 
layer, and which atoms form the lower leaflet, and which 
atoms belong to water molecules. The atom nearest to 
the alcohol molecule then determines in which of the 
three phases a given alcohol molecule is located at any 
given moment. In addition, it is also relevant whether 
an alcohol molecule is hydrogen bonded to some lipid 
molecule, cf. Sec. lIII.Dl 

When all these pieces of information are combined, one 
arrives at data as shown in Fig. 1211 in which we depict 
a few selected ethanol molecules within a DPPC bilayer. 
It is seen that, while the alcohol is inside the bilayer, it is 
hydrogen bonded most of the time. The bonding does not 
persist for the entire duration of a simulation but there 
are short breaks in between. This is in agreement with 
the hydrogen bond lifetime of order 1 ns in Table lTTll while 
ethanol molecules can stay inside the bilayer much longer 
than this. Whenever the hydrogen bond is broken, it 
can either be re-formed shortly afterwards, or the alcohol 
molecule can try to move to some other place. From the 
figure it is seen that an alcohol molecule can move to 
the opposite leaflet of the bilayer but that not all such 
attempts are successful, i.e., the alcohol molecule may 
be reflected back. 

Using the collected information, each alcohol molecule 
is at any given moment in one of five different states (wa- 
ter phase, upper leaflet, upper leaflet hydrogen-bonded, 
lower leaflet and lower leaflet hydrogen-bonded; we will 
discuss the methanol-containing systems a bit further 
down). It is of little interest if an alcohol molecule is 
"scratching" at the surface of the bilayer - rather it is 
important whether the alcohol molecule reaches the part 
of the leaflet where it can form a hydrogen bond. 





successful 


unsuccessful 


System 


number 


time [ps] 


number 


time [ps] 


DPPC-ethanol 


30 


325 


123 


245 


POPC-ethanol 


21 


375 


101 


225 



TABLE V Number of successful and unsuccessful crossing 
events, respectively, within 40 ns of trajectory. In addition, 
the mean time spent in the crossing process is given. 
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FIG. 21 ^-component of the centre-of-mass position of a few tagged ethanol molecule in a DPPC bilayer system. Each figure 
depicts the trajectory of a different molecule. We have chosen three out of the 90 molecules to give a demonstration of the 
possible behaviour. The colours red and blue mark in which leaflet of the bilayer the alcohol is located. If the colour is dark 
red or dark blue, the alcohol is hydrogen bonded to a lipid, otherwise the colour is light red or light blue. If the alcohol is in 
the water phase, the trajectory is shown in yellow. 
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This immediately gives functional definitions for dif- 
ferent kinetic events that can be used for an automatic 
analysis. A successful crossing event from the upper 
to the lower leaflet is for example given by a sequence 
"upper leaflet hydrogen-bonded {upper leaflet} 
{lower leaflet} lower leaflet hydrogen-bonded" where 
the curly braces mean that this step may also be skipped. 
Similarly, an unsuccessful crossing from top to bot- 
tom would be "upper leaflet hydrogen-bonded — > {upper 
leaflet} lower leaflet — > {upper leaflet} upper leaflet 
hydrogen-bonded" . Other criteria arc constructed simi- 
larly. 

Table IVl shows the results of the analysis for crossings 
of ethanol molecules between the two leaflets. A simple 
calculation shows that ethanol molecules are able to move 
from one leaflet to the other on a time scale of 130 ns for 
DPPC and 180 ns for POPC. The number of unsuccessful 
crossing attempts outnumbers the number of successful 
attempts by a factor of 4, thereby demonstrating that 
the hydrophobic tails of the lipid pose a significant bar- 
rier to ethanol not only from the outside of the leaflet 




FIG. 22 View from the top onto part of a POPC bilayer with 
methanol. The lipids are coloured green, methanol is coloured 
blue and a few selected water molecules are shown in red. It 
is easily seen than the methanol is located in a cavity together 
with a few water molecules. 



but also from the inside. For methanol we did not flnd 
any crossing events in our simulations, implying that the 
corresponding time scale must be at least of the order of 
microseconds. 

A closer study of the systems containing methanol is 
hindered by a problem that is not obvious from any of 
the data presented so far. As Fig. 122 shows, methanol 
is virtually never really inside the bilayer, i.e., located 
such that it no longer has direct contact with the bulk 
water phase - in all of our data for DPPC and POPC with 
methanol, we found only a single methanol molecule that 
had actually lost contact with water. 

The above observation is also able to explain why no 
hydrogen bonds are formed between methanol and the 
lipids: Water is simply too energetically favourable a 
binding partner for methanol, or, in other words, the 
surface tension of water is too high for methanol to 
leave the water phase. The concept of an alcohol be- 
ing inside the membrane thus does not apply - topolog- 
ically, methanols are always located outside the mem- 
brane. Rather we need to introduce the concept of a 
methanol being located in a sufficiently deep well. This 
can be quantifled by counting the number of atoms be- 
longing to lipids within a certain distance around some 
particular methanol molecule. This number will be much 
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FIG. 23 Distribution P{t) of the time t for which an alcohol 
molecule stays inside the membrane, left for DPPC, right for 
POPC. Due to limited statistics, we cannot plot P{t) directly 
but are limited to the cumulative probability J P{t')dt' . 
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larger if the methanol is inside such a well. (We use the 
criterion that the number of atoms belonging to lipid 
molecules within 0.6 nm is larger than 50.) 

Using that functional definition, we are able to treat 
ethanol and methanol containing systems on a similar 
footage. While there are no crossing events for methanol, 
another interesting question still arises, namely the dy- 
namics of alcohol exchange between the membrane and 
the water phase. Quantitatively, the interesting quantity 
is the time t between an alcohol molecule entering the 
membrane and its subsequent leaving it again. Our re- 
sults arc shown in Fig. [23 Since there are only of the 
order of 200 (1000) events for ethanol (methanol), the 
statistics is insufficient to compute the probability distri- 
bution P{t). Rather, we present the cumulative probabil- 
ity J P{t')dt', i.e., the probability that an alcohol stays 
inside the membrane no longer than some time t, since 
this quantity can be computed without binning the data 
point. {P{t) follows, in principle, by differentiation of 
the depicted curves.) It is seen from the figure that the 
dynamics is much faster for methanol than for ethanol. 
This comes as no surprise since methanol is not really 
inside the bilayer - it does not need to cross the bilayer 
interface but only needs to deform it (to create a well) . 



IV. DISCUSSION 

This study is, to the best of our knowledge, the first de- 
tailed computational study characterising the behaviour 
of lipid bilayers (POPC and DPPC) under the influence 
of methanol and ethanol. The other existing molecular 
dynamics study of ethanol and POPC (|Feller et al.Ll2002() 
concentrated on the comparisons with an NMR study un- 
der different conditions (at 10 degrees Celsius close to the 
gel phase and using an NVT ensemble). To obtain de- 
tailed information about alcohol-membrane interactions, 
it is thus important to study the biologically important 
fluid phase. 

Let us flrst discuss the area per lipid and bilayer thick- 
ness. The increase in the area per lipid is larger for DPPC 
bilayers (about 7 % for ethanol and 6 % for methanol) 
as compared to POPC systems (5% for ethanol and 
3 % for methanol) . This c ompares well w ith the re- 
cent micropipettc studies of iLv et al.l l)2002() who used 
SOPC vesicles under slightly different conditions (20 vol.- 
% ethanol at room temperature). They observed 9% 
increase in the area per lipid and 8 % decrease in the 
thickness of the bilayer. Here, we obtained a decrease 
of 7 %-10 % (ethanol and DPPC) and 1 %-4 % (ethanol 
POPC) in thickness depending on the definition used, see 
Eq. CEJ. 

As a purely structural effect, it is clear that the mem- 
brane thus becomes more permeable to small molecules 
due to its increased area per lipid. The differences be- 
tween DPPC and POPC are most likely due to the 
slightly longer sn-2 chain of POPC and the double bond 
in it. In addition to the effects captured by the average 



area per lipid, steric constraints seem to make the POPC 
bilayer less susceptible to penetration of small solutes. 

Fur thermore, in a recent study IChanturiva et all 
(|1999() proposed that the penetration of alcohols inside 
the bilayer, and their binding at it, and the resulting de- 
crease in bending rigidity is a feasible pathway for pro- 
moting fusion of cells. Although it is not possible to 
probe this directly by current computer resources, our 
observations support the possibility of such a mechanism. 

It has been suggested that the preferred loca- 
tio n of ethanol close to the mernbrane dehy drates 
it iHolte and GawrischilT99ilKlemmLll99nlll998t> . This 
should show in the water dipole orientation data. Fig. 1181 
The relatively small changes in it, and in the electrostatic 
potential across the membrane, suggest that indirect ef- 
fects, such as receptor blocking, may be more important 
in producing changes in these quantities. 

To characterise thermodynamic properties, we have 
measured the partitioning coefficients. Our estimates of 
K-p for the different alcohol-lipid combinations are given 
in Table IIVI Due to the lack of experimental data, no 
quantitative comparison could be made but qualitatively 
experiments and simulations are in agreement. 

Short-chain alcohols have an amphiphilic character and 
it has been known for long that addition of each new 
CH2 group - adding a CH2 group on methanol gives 
ethanol and so on - has a strong effect on the interactions 
with membran es. This is indicated by the well known 
Traube's rule ijAdamson and Gastl Il997t iTraubd . Il89l|) 
which states that the addition of a new CH2 group leads 
to a decrease in surface tension. In other words, short- 
chain alcohols have a strong effect on membrane proper- 
ties and the effect depends on both the length of the hy- 
drophobic part of the alcohol and on concentra tion. This 
has also been ob served in recent experiments l)Lv et al.l 
I2OO2HLV andTongo. 200^ . 

Our data for methanol and ethanol supports these con- 
clusions. Methanol does not penetrate through the lipid 
tail region which is easily understood by the hydrophobic 
nature of the lipid tails which are repelling methanol as it 
is more polar than ethanol. As a second effect, methanol 
rarely reaches the tail region as each methanol molecule 
moves together with a small cluster of water molecules 
when it is trying to enter the membrane. This means 
that, on one hand, methanol is pulled back into the wa- 
ter phase by this, and, on the other hand, not a single 
small methanol molecule but a significantly larger dressed 
particle, or a small cluster, would need to penetrate the 
membrane. 

The analysis of crossing events, i.e., how often the 
molecules travel through the membrane, showed that 
ethanol is able to penetrate the membrane easily whereas 
for methanol not a single crossing event was observed. 
This confirms the interpretation given above. It is diffi- 
cult to compare these results directly with experiments 
but the possibility of such crossing events has been pro- 
posed on the basis of NMR studies (jHolte and Gawrischl 
Il997j) . The results presented here are, to our knowledge. 
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the first detailed analysis of crossing events. Further ex- 
periments would be needed in order to better characterise 
the situation as the system here is a simple model system 
and the relevance of these results to biological systems, 
in particular yeasts, needs to be better studied. The only 
such a study we were able to find uses NMR and Z. mo- 
bilis (|Schobert et al.lll996fl but direct comparison is not 
possible due to the different experimental setup. 

In the introduction we briefly discussed general 
anaesthesia and membrane-protein interactions induced 
by the addition of anaesthetics, such as small al- 
cohols. ^^^Tliis_was_^bservBd_^n_a,rec^ experi- 
ment l)van den Brink-van der Laan efaH \2004) where 
the potassium channel KcsA was observed to dissoci- 
ate due to the changes in lateral membrane pressure in- 
duced by small alcohols. Here, we have characterised 
simple membrane-alcohol systems. The detailed charac- 
terisation presented here is essential for extensive simula- 
tional studies of membrane-protein-anaesthetic systems. 
From our results it is obvious that the changes in pure 
membranes are subtle but the effects of those changes 
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also sup ported by rec e nt exp eriments using enflurane and 
DPPC ijHauet e t al.|. |2 003tl. S imilar conclusions have 
been drawn bv iTu et alJ |l99^ for the interaction of 
halothane with bilayers. As pointed out by Hauet et al., 
there are various intriguing questions regarding small 
molecules and anaesthesia. These questions are related 
to interactions between membranes and small molecules 
and computer simulations give a direct access to study 
them. 
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